CLAIMS 



We claim: 

1 . A method for processing seismic data to estimate time shift resulting from velocity 
anisotropy in the earth's subsurface, comprising: 

forming a gather of seismic data traces; and 

cross-correlating selected seismic data traces included in said gather within selected 
time windows to estimate the time shift in the seismic data traces included in said gather 
resulting from velocity anisotropy in the earth's subsurface. 

2. The method of claim 1 further comprising adjusting seismic data traces included in 
said gather by the amount of the estimated time shift in said adjusted seismic data trace 
resulting from velocity anisotropy 

3. The method of claim 2 further comprising performing an amplitude variation with 
incidence angle analysis on said adjusted seismic data traces. 

4. The method of claim 2 further comprising performing an amplitude variation with 
azimuth analysis on said adjusted seismic data traces. 

5. The method of claim 2 further comprising: 

determining the incidence angle for seismic data traces included in said gather; and 
applying a least squares analysis process to reflection coefficient, source-receiver 
azimuth angle and incidence angle data for seismic data traces included in said gather to 
calculate the amplitude variation with azimuth and amplitude variation with offset in 
seismic data traces included in said gather. 
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1 6. The method of claim 2 further comprising applying a least squares analysis process 

2 to the time shifts of said seismic data traces to calculate the velocity variation with azimuth 

3 in seismic data traces included in said gather. 

1 7. A method for processing seismic data comprising: 

2 forming a gather of seismic data traces; 

3 performing a surface consistent statics computation on seismic data traces included 

4 in said gather; 

5 cross-correlating successively selected seismic data traces in said gather to estimate 

6 the time shift in seismic data traces included in said gather resulting from azimuthal 

7 velocity anisotropy in the earth's subsurface, applying a least squares analysis process to 

8 the time shifts of said seismic data traces to calculate the velocity variation with azimuth in 

9 seismic data traces included in said gather; 

1 0 utilizing the calculated velocity variations with azimuth to calculate time shifts in 

1 1 seismic data traces included in said gather; 

12 applying said calculated time shifts to said seismic data traces included in said 

13 gather; and 

14 applying a least squares analysis process to reflection coefficient, source-receiver 

1 5 azimuth angle and incidence angle data for seismic data traces included in said gather to 

1 6 calculate the amplitude variation with azimuth and amplitude variation with offset in 

1 7 seismic data traces included in said gather. 

1 8. A method for processing seismic data to estimate time shift resulting from velocity 

2 anisotropy in the earth's subsurface, comprising: 

3 forming a gather of seismic data traces; 

4 forming a pilot trace by combining a selected plurality of said seismic data traces 

5 within a selected time window; and 
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6 cross-correlating a selected seismic data trace included in said gather with said pilot 

7 trace to estimate the time shift in said selected seismic data trace resulting from velocity 

8 anisotropy in the earth's subsurface. 

1 9. A method for processing seismic data to estimate time shift resulting from velocity 

2 anisotropy in the earth's subsurface, comprising: 

3 (a) forming a gather of seismic data traces; 

4 (b) forming a pilot trace by combining a selected plurality of said seismic data 

5 traces within a selected time window; 

6 (c) cross-correlating a selected seismic data trace included in said gather with said 

7 pilot trace to estimate the time shift in said selected seismic data trace resulting from 

8 velocity anisotropy in the earth's subsurface; and 

9 repeating steps (b)and (c) until all seismic data traces within said gather have been 
1 0 cross-correlated with a pilot trace. 

1 1 0. A method for processing seismic data, comprising: 

2 (a) forming a gather of seismic data traces; 

3 (b) forming a pilot trace by combining a selected plurality of said seismic data 

4 traces within a selected time window; 

5 (c) cross-correlating a selected seismic data trace included in said gather with said 

6 pilot trace to estimate the time shift in said selected seismic data trace resulting from 

7 velocity anisotropy in the earth's subsurface; 

8 repeating steps (b) and (c) until all seismic data traces within said gather have been 

9 cross-correlated with a pilot trace; and 

10 adjusting each said selected seismic data trace by the amount of the estimated time 

1 1 shift in each said selected seismic data trace resulting from velocity anisotropy. 

1 11. The method of claim 1 0 further comprising performing an amplitude variation with 

2 incidence angle analysis on said adjusted seismic data traces. 
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1 12. The method of claim 10 further comprising performing an amplitude variation with 

2 azimuth analysis on said adjusted seismic data traces. 

1 13. A method for processing seismic data, comprising: 

2 (a) forming a gather of seismic data traces; 

3 (b) forming a pilot trace by combining a selected plurality of said seismic data 

4 traces within a selected time window; 

5 (c) cross-correlating a selected seismic data trace included in said gather with said 

6 pilot trace to estimate the time shift in said selected seismic data trace resulting from 

7 velocity anisotropy in the earth's subsurface; 

8 repeating steps (b) and (c) until all traces within said gather have been correlated 

9 with a pilot trace; 

1 0 adjusting each said selected seismic data trace by the amount of the estimated time 

1 1 shift in each said selected seismic data trace resulting from velocity anisotropy; 

1 2 determining the incidence angle for each selected seismic data trace; and 

1 3 applying a least squares analysis process to reflection coefficient, source-receiver 

14 azimuth angle and incidence angle data of said seismic data traces to calculate the 

1 5 amplitude variation with azimuth and amplitude variation with offset in seismic data traces 

1 6 included in said gather. 

1 14. A method for processing seismic data, comprising: 

2 (a) forming a gather of seismic data traces; 

3 (b) forming a pilot trace by combining a selected plurality of said seismic data 

4 traces within a selected time window; 

5 (c) cross-correlating a selected seismic data trace included in said gather with said 

6 pilot trace to estimate the time shift in said selected seismic data trace resulting from 

7 velocity anisotropy in the earth's subsurface; 
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8 repeating steps (b) and (c) until all seismic data traces within said gather have been 

9 correlated with a pilot trace; and 

10 applying a least squares analysis process to the time shifts of said seismic data 

1 1 traces to calculate the velocity variation with azimuth in seismic data traces included in 

12 said gather. 

1 15. The method of claim 1 3 further comprising utilizing a least squares analysis to 

2 estimate errors associated with the calculation of amplitude variation in said selected 

3 seismic data traces. 

1 1 6. The method of claim 1 4 further comprising utilizing a least squares analysis to 

2 estimate errors associated with the calculation of time shift variation in said selected 

3 seismic data traces. 

1 1 7. The method of claim 1 4 further comprising utilizing a least squares analysis to 

2 estimate errors associated with the calculation of velocity variation in said selected seismic 

3 data traces. 

1 18. A method for processing seismic data comprising: 

2 (a) forming a gather of seismic data traces; 

3 (b) performing a surface consistent statics computation on said seismic data traces; 

4 (c) cross-correlating successively selected seismic data traces in said gather with a 

5 pilot trace formed by combining a selected plurality of seismic data traces in said gather 

6 within selected time windows to estimate the time shift in said seismic data traces resulting 

7 from azimuthal velocity anisotropy in the earth's subsurface, and applying a least squares 

8 analysis process to the time shifts of said seismic data traces to calculate the velocity 

9 variation with azimuth in seismic data traces included in said gather; 

10 ( d ) utilizing the calculated velocity variations with azimuth to calculate time shifts 

1 1 in seismic data traces included in said gather; 
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12 ( e ) a Pp!ying said calculated time shifts to said seismic data traces included in said 

13 gather; 

14 repeating steps (b) - (e); and 

15 applying a least squares analysis process to reflection coefficient, source-receiver 

16 azimuth angle and incidence angle data of said seismic data traces to calculate the 

1 7 amplitude variation with azimuth and amplitude variation with offset in seismic data traces 

1 8 included in said gather. 

1 19. A method for processing seismic data, comprising: 

2 (a) forming a gather of seismic data traces; 

3 (b) selecting a time window of said gather; 

4 (c) selecting a plurality of said traces; 

5 (d) summing said plurality of traces within said time widow; 

6 (e) selecting a trace from said gather and cross-correlating said selected trace with 

7 said summed plurality of traces; 

8 (f) determining a time displacement of said selected trace which results in a 

9 substantially maximum correlation of said selected trace with said summed plurality of 

1 0 traces within said selected time window; 

1 1 (g) repeating steps (c) - (f) for successive traces and selected summed pluralities of 

12 traces until a time displacement has been determined for each trace within said gather 

1 3 which results in a substantially maximum correlation of each trace of said gather with a 

14 summed plurality of traces within said time window; and 

15 (h) repeating steps (b) - (g) for successive time windows within said gather. 

1 20. A method for processing seismic data to compensate for time shifts in said data 

2 resulting from azimuthal anisotropy, comprising: 

3 (a) forming a gather of seismic data traces; 

4 (b) selecting a time window of said gather; 

5 (c) selecting a plurality of said traces; 
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6 (d) summing said plurality of traces within said time widow; 

7 (e) selecting a trace from said gather and cross-correlating said selected trace with 

8 said summed plurality of traces; 

9 (f) determining a time displacement of said selected trace which results in a 

10 substantially maximum correlation of said selected trace with said summed plurality of 

1 1 traces within said selected time window; 

12 (g) repeating steps (c) - (f) for successive traces and selected summed pluralities of 

13 traces until a time displacement has been determined for each trace within said gather 

14 which results in a substantially maximum correlation of each trace of said gather with a 

1 5 summed plurality of traces within said time window; 

16 (h) repeating steps (b) - (g) for successive selected time windows of said gather; 

17 (i) applying the determined time displacement for each trace within each said 

1 8 window to the trace at the center point of each said time window; and 

19 (j) interpolating time shifts for the remaining portion of each data trace between the 

20 window center points. 

1 21. The method of claim 20 further comprising: 

2 determining the reflection coefficient and azimuthal angle for each trace of the 

3 gather; 

4 calculating the value of the incidence angle for each trace of the gather; and 

5 utilizing a least squares method to compute reflection coefficients as a function of 

6 azimuthal angle and incident angle for the seismic traces of the gather. 

1 22. The method of claim 20 further comprising utilizing a linear least squares method 

2 to calculate the semi-major and semi-minor axis, and the angle between a predefined zero 

3 azimuth direction and the semi-major axis direction of the an ellipse representing the 

4 azimuthal variation in velocity. 
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1 23. The method of claim 22 further comprising: 

2 formulating a matrix for a least squares analysis for estimating the errors associated 

3 with the calculation of the time shift variation with azimuth; 

4 computing an error estimate by taking the square root of the diagonals of said 

5 matrix; and 

6 comparing the estimated error to the computation of the azimuthal velocity 

7 variation to confirm whether or not the azimuthal velocity variation has an acquisition 

8 footprint. 

1 24. A digital computer programmed to utilize seismic data traces obtained over a region 

2 of the earth's subsurface to perform a process comprising the steps of: 

3 forming a gather of seismic data traces; 

4 forming a pilot trace by combining a selected plurality of said seismic data traces 

5 within a selected time window; and 

6 cross-correlating a selected seismic data trace included in said gather with said pilot 

7 trace to estimate the time shift in said selected seismic data trace resulting from velocity 

8 anisotropy in the earth's subsurface. 

1 25. A device which is readable by a digital computer having instructions defining the 

2 following process and instructions to the computer to perform said process: 

3 forming a gather of seismic data traces; 

4 forming a pilot trace by combining a selected plurality of said seismic data traces 

5 within a selected time window; and 

6 cross-correlating a selected seismic data trace included in said gather with said pilot 

7 trace to estimate the time shift in said selected seismic data trace resulting from velocity 

8 anisotropy in the earth's subsurface. 

1 26. The device of claim 25 wherein said device is selected from the group consisting of 

2 a magnetic tape a magnetic disk, an optical disk and a CD-ROM. 
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